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The development of primordial inhomogeneities into the non-linear regime and 
the formation of the first astrophysical objects within dark matter halos mark 
the transition from a simple, neutral, cooling universe - described by just a few 
parameters - to a messy ionized one — the realm of radiative, hydrodynamic, and 
star formation processes. The recent measurement by the WMAP satellite of a 
large optical depth to electron scattering implies that this transition must have 
begun very early, and that the universe was reionized at redshift z- lon = 17 ± 5. It 
is an early generation of extremely metal-poor massive stars and/or 'seed' accreting 
black holes in subgalactic halos that may have generated the ultraviolet radiation 
and mechanical energy that reheated and reionized most of the hydrogen in the 
cosmos. The detailed thermal, ionization, and chemical enrichment history of the 
universe during the crucial formative stages around z = 10 — 20 depends on the 
power-spectrum of density fluctuations on small scales, the stellar initial mass 
function and star formation efficiency, a complex network of poorly understood 
'feedback' mechanisms, and remains one of the crucial missing links in galaxy 
formation and evolution studies. 



1. Introduction 

The last decade has witnessed great advances in our understanding of the 
high redshift universe. The pace of observational cosmology and extra- 
galactic astronomy has never been faster, and progress has been equally 
significant on the theoretical side. The key idea of currently popular cos- 
mological scenarios, that primordial density fluctuations grow by gravi- 
tational instability driven by cold, collisionless dark matter (CDM), has 
been elaborated upon and explored in detail through large-scale numerical 
simulations on supercomputers, leading to a hierarchical ('bottom-up') sce- 
nario of structure formation. In this model, the first objects to form are on 
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subgalactic scales, and merge to make progressively bigger structures ('hier- 
archical clustering'). Ordinary matter in the universe follows the dynamics 
dictated by the dark matter until radiative, hydrodynamic, and star forma- 
tion processes take over. Perhaps the most remarkable success of this theory 
has been the prediction of anisotropics in the temperature of the cosmic 
microwave background (CMB) radiation at about the level subsequently 
measured by the COBE satellite and most recently by the BOOMERANG, 
MAXIMA, DASI, CBI, Archeops, and WMAP experiments. 

In spite of some significant achievements in our understanding of the 
formation of cosmic structures, there are still many challenges facing hi- 
erarchical clustering theories, and many fundamental questions remain, at 
best, only partially answered. While quite successful in matching the ob- 
served large-scale density distribution (like, e.g., the properties of galaxy 
clusters, galaxy clustering, and the statistics of the Lyman-a forest), CDM 
simulations appear to produce halos that are too centrally concentrated 
compared to the mass distribution inferred from the rotation curves of 
(dark matter-dominated) dwarf galaxies, and to predict too many dark mat- 
ter subhalos compared to the number of dwarf satellites observed within 
the Local Group. 38 ' 53 ' 37 ' 27 Another perceived problem (possibly connected 
with the 'missing satellites' 8 ) is our inability to predict when, how, and 
to what temperature the universe was reheated and reionized, i.e. to un- 
derstand the initial conditions of the galaxy formation process. While N- 
body+hydrodynamical simulations have convincingly shown that the inter- 
galactic medium (IGM) the main repository of baryons at high rcdshift 
- is expected to fragment into structures at early times in CDM cosmogo- 
nies, the same simulations arc much less able to predict the efficiency with 
which the first gravitationally collapsed objects lit up the universe at the 
end of the 'dark ages'. The crucial processes of star formation, preheating 
and feedback (e.g. the effect of the heat input from the first generation of 
sources on later ones), and assembly of massive black holes in the nuclei 
of galaxies are poorly understood. 30 We know that at least some galaxies 
and quasars were already shining when the universe was less than 10 9 yr 
old. But when did the first luminous objects form, what was their nature, 
and what impact did they have on their environment and on the forma- 
tion of more massive galaxies? While the excess H i absorption measured 
in the spectra of z <~ 6 quasars in the Sloan Digital Sky Survey (SDSS) 
has been interpreted as the signature of the trailing edge of the cosmic 
reionization epoch 3 ' 17 ' 15 , the recent detection by the Wilkinson Microwave 
Anisotropy Probe ( WMAP) of a large optical depth to Thomson scatter- 
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ing, r e = 0.17 ± 0.04 suggests that the universe was reionized at higher 
redshifts, z- lon = 17 ± 5. 28 ' 51 This is of course an indication of significant 
star- format ion activity at very early times. 

In this talk I will summarize some recent developments in our under- 
standing of the dawn of galaxies and the impact that some of the earliest 
cosmic structure may have had on the baryonic universe. 

2. The dark ages 

In the era of precision cosmology we know that, at a redshift Zd cc = 1088±1, 
exactly idee = (372 ± 14) x 10 3 years after the big bang, the universe be- 
came optically thin to Thomson scattering 51 , and entered a 'dark age'. 42 
At this epoch the electron fraction dropped below 13% (Figure 1), and 
the primordial radiation cooled below 3000 K, shifting first into the in- 
frared and then into the radio. We understand the microphysics of the 
post-recombination universe well. The fractional ionization froze out to 
the value ~ 10 _4 ' 8 f2M / '(M)&): these residual electrons were enough to keep 
the matter in thermal equilibrium with the radiation via Compton scat- 
tering until a thermalization redshift z t ~ 800(f2{,/i 2 ) 2 / 5 ~ 150, i.e. well 
after the universe became transparent. 40 Thereafter, the matter temper- 
ature decreased as (1 + z) 2 due to adiabatic expansion (Figure 2) until 
primordial inhomogeneities in the density field evolved into the non-linear 




Figure 1. Helium and hydrogen recombination for the WMAP parameters 
(U M ,n A ,n b ,h) = (0.29, 0.71,0.045, 0. 7). 51 The step at earlier times in the left panel 
is due to the recombination of He III into He II . We used the code RECFAST 47 to 
compute the electron fraction. 
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regime. The minimum mass scale for the gravitational aggregation of cold 
dark matter particles is negligibly small. One of the most popular CDM 
candidates is the neutralino: in neutralino CDM, collisional damping and 
free streaming smear out all power of primordial density inhomogeneities 
only below ~ 10~ 7 Mq. 26 Baryons, however, respond to pressure gradients 
and do not fall into dark matter clumps below the cosmological Jeans mass 
(in linear theory this is the minimum mass-scale of a perturbation where 
gravity overcomes pressure), 

Mj = *f (Ig^) ^ - 2.5 x 10^- M (aT/^Q M V 2 . (1) 

Here a = (1 + is the scale factor, p the total mass density including 
dark matter, /i the mean molecular weight, and T the gas temperature. The 
evolution of Mj is shown in Figure 2. In the post-recombination universe, 
the baryon-electron gas is thermally coupled to the CMB, T oc a" 1 , and 
the Jeans mass is independent of redshift and comparable to the mass of 
globular clusters, Mj w 10 6 M©. For z < Zt, the temperature of the 
baryons drops as T oc a" 2 , and the Jeans mass decreases with time, Mj oc 
a~ 3 / 2 . This trend is reversed by the reheating of the IGM. The energy 
released by the first collapsed objects drives the Jeans mass up to galaxy 
scales (Figure 2): previously growing density perturbations decay as their 




Figure 2. Left: Evolution of the radiation (long-dashed line, labeled CMB) and gas 
(solid line, labeled IGM) temperatures after recombination. The universe is assumed to 
be reionized by ultraviolet radiation at z ~ 20. The short-dashed line is the extrapolated 
gas temperature in the absence of any reheating mechanism. Right: Cosmological (gas 
+ dark matter) Jeans (solid line) and filtering (dot-dashed line) mass. 
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mass drops below the new Jeans mass. In particular, photoionization by the 
ultraviolet radiation from the first stars and quasars would heat the IGM 
to temperatures of w 10 4 K (corresponding to a Jeans mass Mj < 10 10 M Q 
at z ~ 20), suppressing gas infall into low mass halos and preventing new 
(dwarf) galaxies from forming. 

3. Linear theory 

When the Jeans mass itself varies with time, linear gas fluctuations tend to 
be smoothed on a (filtering) scale that depends on the full thermal history of 
the gas instead of the instantaneous value of the sound speed. 23 From linear 
perturbation analysis, and for a flat universe at high redshift, the growth 
of density fluctuations in the gas is suppressed for comoving wavenumbers 
k > kp, where the filtering scale ftp is related to the Jeans wavenumber kj 



Here kj = {a/ c s )^J AnGp, and c s is the sound speed. This expression for 
kp accounts for an arbitrary thermal evolution of the IGM through kj(a). 
Corresponding to the critical wavenumber kp there is a critical (filtering) 
mass Mp, defined as the mass enclosed in the sphere with comoving radius 
equal to kp, 



The Jeans mass Mj is defined analogously in terms of kj. It is the fil- 
tering mass that is central to calculations of the effects of reheating and 
reionization on galaxy formation. The filtering mass for a toy model with 
early photoionization is shown in Figure 2: after reheating, the filtering 
scale is actually smaller than the Jeans scale. Numerical simulations of 
cosmological reionization confirm that the characteristic suppression mass 
is typically lower than the linear-theory Jeans mass. 22 

4. The emergence of cosmic structure 

As mentioned in the introduction, some shortcomings on galactic and sub- 
galactic scales of the currently favored model of hierarchical galaxy forma- 
tion in a universe dominated by CDM have recently appeared. The signif- 
icance of these discrepancies is still debated, and 'gastrophysicaP solutions 
involving feedback mechanisms may offer a possible way out. Other models 
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Alp = (47r/3)p(27ra/fcF) 3 . 
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have attempted to solve the apparent small-scale problems of CDM at a 
more fundamental level, i.e. by reducing small-scale power. Although the 
'standard' ACDM model for structure formation assumes a scale-invariant 
initial power spectrum of density fluctuations, P(k) cx k n with n = 1, the 
recent WMAP data favor (but don't require) a slowly varying spectral in- 
dex, dn/d\uk = — 0.031±g;gi|, i.e. a model in which the spectral index 
varies as a function of wavenumber k. 51 This running spectral index model 
predicts a significanly lower amplitude of fluctuations on small scales than 
standard ACDM. The suppression of small-scale power has the advantage of 
reducing the amount of substructure in galactic halos and makes small halos 
form later (when the universe was less dense) hence less concentrated, 38 ' 57 
relieving some of the problems of ACDM. But it makes early reionization 
a challenge. 

Figure 3 shows the linearly extrapolated (to 2 = 0) variance of the 



Figure 3. The variance of the matter-density field vs. mass M, for different power spec- 
tra. All models assume a 'concordance' cosmology with parameters (Qm > ^A> ^f>> h) = 
(0.29,0.71,0.045,0.7). Solid curve: standard ACDM with no tilt, cluster normalized. 
Dotted curve: AWDM with a particle mass mx = 2keV, cluster normalized, no tilt. 
Dashed curve: tilted WMAP model, WMAP data only. Dash-dotted curve: tilted 
WMAP model, including 2dFGRS and Lyman-a data. Dash-triple dotted curve: run- 
ning spectral index WMAP model, including 2dFGRS and Lyman-o data. Here n refers 
to the spectral index at k = 0.05 Mpc -1 . The horizontal line at the top of the figure 
shows the value of the extrapolated collapse overdensity S c (z) at z = 20. 
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mass-density field smoothed on a scale of comoving radius R, 

i r°° 

a 2 M = ((SM/M) 2 ) = -^ dkk 2 P{k)T 2 (k)W 2 (kR), (4) 
ZlT Jo 

for different power spectra. Here M — HqUmR 3 /2G is the mass inside R, 
T{k) is the transfer function for the matter density field (which accounts for 
all modifications of the primordial power-law spectrum due to the effects 
of pressure and dissipative processes), and W(kR) is the Fourier transform 
of the spherical top-hat window function, W(x) = (3/x 2 )(sinx/x — cosx). 
The value of the rms mass fluctuation in a 8 Mpc sphere, as = cr(z — 
0, R = 8/i _1 Mpc), has been fixed for the n = 1 models to cr$ = 0.74, 
consistent with recent normalization by the z — X-ray cluster abundance 
constraint. 43 

In the CDM paradigm, structure formation proceeds 'bottom-up', i.e., 
the smallest objects collapse first, and subsequently merge together to 
form larger objects. It then follows that the loss of small-scale power 
modifies structure formation most severely at the highest redshifts, sig- 
nificantly reducing the number of self-gravitating objects then. This, of 
course, will make it more difficult to reionize the universe early enough. It 
has been argued, for example, that one popular modification of the CDM 
paradigm, warm dark matter (WDM), has so little structure at high red- 
shift that it is unable to explain the WMAP observations of an early epoch 
of reionization. 51,2 And yet the WMAP running-index model may suffer 
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Figure 4. Mass fraction in all collapsed halos above the filtering mass Mp as a function 
of redshift, for different power spectra. Curves are the same as in Figure 3. Left panel: 
filtering mass Mp has been computed in the absence of reionization. Right panel: Mp 
computed assuming the universe is reionized by ultraviolet radiation at z ~ 20. 
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from a similar problem. 49 A look at Figure 3 shows that 10 6 M halos will 
collapse at z — 20 from 2.9 a fluctuations in a tilted ACDM model with 
n = 0.99 and cr 8 = 0.9, from 4.6 a fluctuations in a running-index model, 
and from 5.7 c fluctuations in a WDM cosmology. The problem is that 
scenarios with increasingly rarer halos at early times require even more ex- 
treme assumptions (i.e. higher star formation efficiencies and UV photon 
production rates) in order to be able to reionize the universe by z ~ 17 
as favored by WMAP. 49 ' 56 ' 24 ' 11 ' 9 Figure 4 depicts the mass fraction in all 
collapsed halos with masses above the filtering mass for a case without 
reionization and one with reionization occuring at z ~ 20. At early epochs 
this quantity appears to vary by orders of magnitude in different models! 

5. The epoch of reionization 

Since hierarchical clustering theories provide a well-defined framework in 
which the history of baryonic material can be tracked through cosmic time, 
probing the reionization epoch may then help constrain competing mod- 
els for the formation of cosmic structures. Quite apart from uncertainties 
in the primordial power spectrum on small scales, however, it is the as- 
trophysics of baryons that makes us unable to predict when reionization 
actually occurred. Consider the following illustrative example: 

Hydrogen photoionization requires more than one photon above 13.6 
eV per hydrogen atom: of order t/t rcc ~ 10 (where t rcc is the volume- 
averaged hydrogen recombination timescale) extra photons appear to be 
needed to keep the gas in overdense regions and filaments ionized against 
radiative recombinations. 21 ' 33 . A 'typical' stellar population produces dur- 
ing its lifetime about 4000 Lyman continuum (ionizing) photons per stellar 
proton. A fraction / ~ 0.25% of cosmic baryons must then condense into 
stars to supply the requisite ultraviolet flux. This estimate assumes a stan- 
dard (Salpetcr) initial mass function (IMF), which determines the relative 
abundances of hot, high mass stars versus cold, low mass ones. 

The very first generation of stars ('Population III') must have formed, 
however, out of unmagnetized metal-free gas: numerical simulations of the 
fragmentation of pure H and He molecular clouds 6 ' 1 have shown that these 
characteristics likely led to a 'top-heavy' IMF biased towards very mas- 
sive stars (VMSs, i.e. stars a few hundred times more massive than the 
Sun), quite different from the present-day Galactic case. Metal- free VMSs 
emit about 10 5 Lyman continuum photons per stellar baryon 7 , approxi- 
mately 25 times more than a standard stellar population. A corresponding 
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smaller fraction of cosmic baryons would have to collapse then into VMSs 
to reionize the universe, / ~ 10~ 4 . There are of course further complica- 
tions. Since, at zero metallicity, mass loss through radiatively-driven stellar 
winds is expected to be negligible 29 , Population III stars may actually die 
losing only a small fraction of their mass. If they retain their large mass 
until death, VMSs with masses 100 < m < 250 M will encounter the 
electron-positron pair instability and disappear in a giant nuclear-powered 
explosion 18 , leaving no compact remnants and polluting the universe with 
the first heavy elements. In still heavier stars, however, oxygen and silicon 
burning is unable to drive an explosion, and complete collapse to a black 
hole will occur instead. 5 Thin disk accretion onto a Schwarzschild black 
hole releases about 50 MeV per baryon. The conversion of a trace amount 
of the total baryonic mass into early black holes, / ~3x 10~ 6 , would then 
suffice to reionize the universe. 

6. Preheating and galaxy formation 

Even if the IMF at early times were known, we still would remain uncer- 
tain about the fraction of cold gas that gets retained in protogalaxies after 
the formation of the first stars (this quantity affects the global efficiency 
of star formation at these epochs) and whether - in addition to ultraviolet 
radiation - an early input of mechanical energy may also play a role in de- 
termining the thermal and ionization state of the IGM on large scales. The 
same massive stars that emit ultraviolet light also explode as supernovae 
(SNe), returning most of the metals to the interstellar medium of pregalac- 
tic systems and injecting about 10 51 ergs per event in kinetic energy. A 
complex network of feedback mechanisms is likely at work in these sys- 
tems, as the gas in shallow potential is more easily blown away, 13 thereby 
quenching star formation. Furthermore, as the blastwaves produced by 
supernova explosions - and possibly also by winds from 'miniquasars' - 
sweep the surrounding intergalactic gas, they may inhibit the formation 
of nearby low-mass galaxies due to 'baryonic stripping' 45 , and drive vast 
portions of the IGM to a significantly higher temperature than expected 
from photoionization, 54,31,32,52 ' 10 so as to 'choke off' the collapse of fur- 
ther galaxy-scale systems. Note that this type of global feedback is fun- 
damentally different from the 'in situ' heat deposition commonly adopted 
in galaxy formation models, in which hot gas is produced by supernovae 
within the parent galaxy. We refer here to this global early energy input 
as 'preheating'. 4 Note that a large scale feedback mechanism may also be 
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operating in the intraclustcr medium: studies of X-ray emitting gas in clus- 
ters show evidence for some form of non-gravitational entropy input 41 . 
The energy required there is at a level of ~ 1 keV per particle, and must be 
injected either in a more localized fashion or at late epochs in order not to 
violate observational constraints on the temperature of the Lyman-a forest 
at z ~ 3. The thermal and ionization history of a preheated universe may 
be very different from one where hydrogen is photoionized. The gas will be 
heated up to a higher adiabat, and collisions with hot electrons will be the 
dominant ionization mechanism. The higher energies associated with pre- 
heating may doubly ionize helium at high-z, well before the 'quasar epoch' 
at z ~ 3. Galaxy formation and evolution will also be different, as pre- 
heating will drive the filtering mass above 10 10 — 10 11 M Q and will tend to 
flatten the faint-end slope of the present-epoch galaxy luminosity function, 
in excellent agreement with the data and without the need for SN feedback 
at late times. 4 

It is interesting to set some general constraints on the early star- 
formation episode and stellar populations that may be responsible for an 
early preheating of the IGM at the levels consistent with the temperature 
of intergalactic gas inferred at z ~ 3. Let us characterize the energy input 
due to preheating by the energy per baryon, E p , deposited in the IGM at 
redshift z p . We examine a homogenous energy deposition since the filling 
factor of pregalactic outflows is expected to be large. 32 ' 19 Let fi* be the 
mass density of stars formed at z p in units of the critical density, i?sN the 
mechanical energy injected per SN event, and f w the fraction of that energy 
that is eventually deposited into the IGM. Denoting with r\ the number of 
SN explosions per mass of stars formed, one can write 



where m p is the proton mass. For a Salpeter IMF between 0.1 and 100 M , 
the number of Type II SN explosions per mass of stars formed is 77 = 0.0074 
Mq 1 , assuming all stars above 8 M© result in SNc II. Numerical simulations 
of the dynamics of SN-driven bubbles from subgalactic halos have shown 
that up to 40% of the available SN mechanical luminosity can be converted 
into kinetic energy of the blown away material, f w w 0.4, the remainder 
being radiated away. 36 With E$n = 1-2 x 10 51 ergs, equation (5) then 
implies 
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SN-driven pregalactic outflows efficiently carry metals into intergalactic 
space. 32 For a normal IMF, the total amount of metals expelled in winds 
and final ejecta (in SNe or planetary nebulae) is about 1% of the input 
mass. Assuming a large fraction, fz = 0.5, of the metal-rich SN ejecta es- 
cape the shallow potential wells of subgalactic systems, the star-formation 
episode responsible for early preheating will enrich the IGM to a mean level 

= °- 01 "* /z = 0.014 Z Q (E p /0.l keV). (7) 
Mb 

The weak C IV absorption lines observed in the Lyman-a forest at z = 
3 — 3.5 imply a minimum universal mctallicity relative to solar in the range 
[-3.2] to [-2.5]. 50 . Preheating energies in excess of 0.1 keV appear then 
to require values of fi» and (Z) that are too high, comparable to the total 
mass fraction in stars seen today 20 and in excess of the enrichment of the 
IGM inferred at intermediate redshifts, respectively. 

The astrophysics of first light may not be as simple, however. The metal 
constraint assumes that metals escaping from protogalaxics are evenly 
mixed into the IGM and the Lyman-a clouds. 54 Inefficient mixing could 
instead produce a large variance in intergalactic mctallicities. The metal 
abundances of the Lyman-a clouds may underestimate the average metallic- 
ity of the IGM if there existed a significant warm-hot gas phase component 
with a higher level of enrichment, as detected for example in O VI. 48 Today, 
the metallicity of the IGM may be closer to ~ 1/3 of solar if the metal pro- 
ductivity of galaxies within clusters is to be taken as representative of the 
universe as a whole. 44 Uncertainties in the early IMF make other preheat- 
ing scenarios possible and perhaps even more likely. Population III stars 
with main-sequence masses of approximately 140 — 260 M Q will encounter 
the electron-positron pair instability and be completely disrupted by a gi- 
ant nuclear-powered explosion. 25 A fiducial 200 M Q Population III star will 
explode with a kinetic energy at infinity of E$n = 4 x 10 52 ergs, injecting 
about 90 Mq of metals. For a very 'top-heavy' IMF with rj = 0.005 M" 1 , 
equation (5) now yields 

= 0.001 (Sp/O.lkeV), (8) 

in 

and a mean IGM metallicity 

= °' 45 "* fz = 0.02 Z Q (E p /0.1 keV) (9) 
lib 

(in both expressions above we have assumed f w = fz = This scenario 
can yield large preheating energies by converting only a small fraction of 
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the comic baryons into Population III stars, but tends to produce too many 
metals for E p > 0.1 keV. The metallicity constraint, of course, does not 
bound preheating from winds produced by an early, numerous population 
of faint 'miniquasars'. a Accretion onto black holes releases 50 MeV per 
baryon, and if a fraction f w of this energy is used to drive an outflow and is 
ultimately deposited into the IGM, the accretion of a trace amount of the 
total baryonic mass onto early black holes, 

may then suffice to preheat the whole universe. Note that this value is about 
50 f w times smaller than the density parameter of the supermassive variety 
found today in the nuclei of most nearby galaxies, I^smbh ~ 2 x 10~ 6 /i -1 . 35 



7. Conclusions 

The above discussion should make it clear that, despite much recent 
progress in our understanding of the formation of early cosmic structure 
and the high-redshift universe, the astrophysics of first light remains one 
of the missing links in galaxy formation and evolution studies. We are 
left very uncertain about the whole era from 10 8 to 10 9 yr - the epoch 
of the first galaxies, stars, supernovae, and massive black holes. Some of 
the issues discussed above are likely to remain a topic of lively controversy 
until the launch of the James Webb Space Telescope (JWST), ideally suited 
to image the earliest generation of stars in the universe. If the first mas- 
sive black holes form in pregalactic systems at very high redshifts, they 
will be incorporated through a series of mergers into larger and larger ha- 
los, sink to the center owing to dynamical friction, accrete a fraction of 
the gas in the merger remnant to become supermassive, and form binary 
systems. 55 Their coalescence would be signalled by the emission of low- 
frequency gravitational waves detectable by the planned Laser Interferom- 
eter Space Antenna (LISA). An alternative way to probe the end of the dark 
ages and discriminate between different reionization histories is through 21 
cm tomography. 34 Prior to the epoch of full reionization, 21 cm spectral fea- 
tures will display angular structure as well as structure in redshift space due 

a Because the number density of bright quasi-stellar objects at z > 3 is low 16 , the thermal 
and kinetic energy they expel into intergalactic space must be very large to have a global 
effect, i.e. for their blastwaves to fill and preheat the universe as a whole. The energy 
density needed for rare, luminous quasars to shock-heat the entire IGM would in this 
case violate the COBE limit on y-distortion. 54 
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to inhomogcneities in the gas density field, hydrogen ionized fraction, and 
spin temperature. Radio maps will show a patchwork (both in angle and 
in frequency) of emission signals from H i zones modulated by H n regions 
where no signal is detectable against the CMB. 12 The search at 21 cm for 
the epoch of first light has become one of the main science drivers of the 
LOw Frequency ARray (LOFAR). While remaining an extremely challeng- 
ing project due to foreground contamination from unresolved extragalactic 
radio sources 14 and free-free emission from the same halos that reionize 
the universe 39 , the detection and imaging of large-scale structure prior to 
reionization breakthrough remains a tantalizing possibility within range of 
the next generation of radio arrays. 

Acknowledgments 

We have benefited from many discussions with all our collaborators on the 
topics described here. PM acknowledges the hospitality of the Carnegie 
Observatories where part of this review was written. Support for this work 
was provided by NSF grant AST-0205738 and by NASA grant NAG5-11513. 

References 

1. Abel, T., Bryan, G., & Norman, M. 2000, ApJ, 540, 39 

2. Barkana, R., Haiman, Z., & Ostriker, J. P. 2001, ApJ, 558, 482 

3. Becker, R. H., et al.2001, AJ, 122, 2850 

4. Benson, A. J., & Madau, P. 2003, MNRAS, submitted (astro-ph/0303121) 

5. Bond, J. R., Arnett, W. D., & Carr, B. J. 1984, ApJ, 280, 825 

6. Bromm, V., Coppi, P. S., & Larson, R. B. 2002, ApJ, 564, 23 

7. Bromm, V., Kudritzki, R. P., & Loeb, A. 2001, ApJ, 552, 464 

8. Bullock, J. S., Kravtsov, A. V., & Weinberg, D. H. 2001, ApJ, 548, 33 

9. Cen, R. 2003, ApJ, submitted (astro-ph/0303236) 

10. Cen, R., & Bryan, G. L. 2001, ApJ, 546, L81 

11. Ciardi, B., Ferrara, A., & White, S. D. M. 2003, MNRAS, submitted (astro- 
ph/0302451) 

12. Ciardi, B., & Madau, P. 2003, ApJ, submitted (astro-ph/0303249) 

13. Dekel, A., & Silk, J. 1986, ApJ, 303, 39 

14. Di Matteo, T., Perna, R., Abel, T., & Rees, M. J. 2002, ApJ, 564, 576 

15. Djorgovski, S. G., Castro, S. M., Stern, D., & Mahabal, A. A. 2001, ApJ, 
560, L5 

16. Fan, X., et al.2001, AJ, 121, 54 

17. Fan, X., et al.2002, AJ, 123, 1247 

18. Fryer, C. L., Woosley, S. E., & Heger, A. 2001, ApJ, 550, 372 

19. Furlanetto S., & Loeb A. 2003, ApJ, in press (astro-ph/0211496) 

20. Glazebrook, K., et al.2003, ApJ, in press (astro-ph/0301005) 



14 



21. Gnedin, N. Y. 2000, ApJ, 535, 530 

22. Gnedin, N. Y. 2000, ApJ, 542, 535 

23. Gnedin, N. Y., & Hui, L. 1998, MNRAS, 296, 44 

24. Haiman, Z., & Holder, G. P. 2003, ApJ, submitted (astro-ph/0302403) 

25. Heger, A., & Woosley, S. E. 2002, ApJ, 567, 532 

26. Hofmann, S., Schwarz, D. J., & Stocker, H. 2001, PhRvD, 64, 083507 

27. Klypin, A., Kravtsov, A. V., Valenzuela, O., & Prada, F. 1999, ApJ, 522, 82 

28. Kogut, A., et al.2003, ApJ, submitted (astro-ph/0302213) 

29. Kudritzki, R. P. 2000, in The First Stars, ed. A. Weiss, T. Abel, & V. Hill 
(Heidelberg: Springer), 127 

30. Locb, A., & Barkana, R. 2001, ARA&A, 39, 19 

31. Madau, P. 2000, RSPTA, 358, 2021 

32. Madau, P., Ferrara, A., & Rees, M. J. 2001, ApJ, 555, 92 

33. Madau, P., Haardt, F., & Rees, M. J. 1999, ApJ, 514, 648 

34. Madau, P., Meiksin, A., & Rees, M. J. 1997, ApJ, 475, 492 

35. Merritt D., & Ferrarese, L. 2001, ApJ, 547, 140 

36. Mori, M., Ferrara, A., & Madau, P. 2002, ApJ, 571, 40 

37. Moore, B., et al.1999, ApJ, 524, L19 

38. Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 

39. Oh, S. P., & Mack, K. J. 2003, MNRAS, submitted (astro-ph/0302099) 

40. Peebles, P. J. E. 1993, Principles of Physical Cosmology (Princeton: Prince- 
ton University Press) 

41. Ponman, T. J., Cannon, D. B., & Navarro, J. F. 1999, Nature, 397, 135 

42. Rees, M. J. 2000, RSPTA, 358, 1989 

43. Reiprich, T. H., & Bohringer, H. 2002, ApJ, 567, 716 

44. Renzini, A. 1997, ApJ, 488, 35 

45. Scannapicco, E., Ferrara, A., & Madau, P. 2002, ApJ, 574, 590 

46. Schneider, R., Ferrara, A., Natarajan, P., & Omukai, K. 2002, ApJ, 571, 30 

47. Seager, S., Sasselov, D. D., & Scott, D. 1999, ApJ, 523, LI 

48. Simcoe, R. A., Sargent, W. L. W., & Rauch, M. 2002, ApJ, 578, 737 

49. Somerville, R. S., Bullock, J. S., & Livio, M. 2003, ApJ, submitted (astro- 
ph/0303481) 

50. Songaila, A. 1997, ApJ, 490, LI 

51. Spcrgel, D.N., et al.2003, ApJ, submitted (astro-ph/0302207) 

52. Theuns, T., Mo, H.J., & Schaye, J. 2001, MNRAS, 321, 450 

53. van den Bosch, F. C, Robertson, B. E., Dalcanton, J., & de Block, W. J. G. 
2000, AJ, 119, 1579 

54. Voit, G. M. 1996, ApJ, 465, 548 

55. Volonteri, M., Haardt, F., & Madau, P. 2003, ApJ, 582, 559 

56. Wyithe, J. S. B., & Loeb, A. 2003, ApJ, submitted (astro-ph/0302297) 

57. Zentner, A. R., & Bullock, J. S. 2002, PhRvD, 66, 43003 



